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ABSTRACT 

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on mea¬ 
surements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar 
masses and star formation rates of six high-redshift (2.8 < 2 ; < 5.7) dusty, star-forming galaxies 
(DSEGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first 
discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic 
redshifts and robust lens models derived from ALMA observations. We have conducted follow-up ob¬ 
servations, obtaining multi-wavelength imaging data, using HST^ Spitzer^ Herschel and the Atacama 
Pathfinder Experiment (APEX). We use the high-resolution HST/WFC3 images to disentangle the 
background source from the foreground lens in Spitzer/IRAC data. The detections and upper limits 
provide important constraints on the spectral energy distributions (SEDs) for these DSEGs, yielding 
stellar masses, IR luminosities, and star formation rates (SERs). The SED fits of six SPT sources 
show that the intrinsic stellar masses span a range more than one order of magnitude with a median 
value ^ 5 xIO^^Mq. The intrinsic IR luminosities range from dxIO^^LQ to 4x10 ^^Lq. They all have 
prodigious intrinsic star formation rates of 510 to 4800 Gompared to the star-forming main 

sequence (MS), these six DSEGs have specific SERs that all lie above the MS, including two galaxies 
that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing 
strong starburst events which may be driven by major mergers. 

Subject headings: galaxies: high-redshift 


1. INTRODUCTION 

Studies of the Gosmic Infrared Background have shown 
that about half of the energy produced in cosmic his¬ 
tory comes from distant dust-enshrouded star-forming 
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galaxies (jPole et al.l l2QQ6[ ). Ever since the first extra- 
galac t ic snbmillimeter-wavelength surveys (jSmail et al.l 
I 1997 I : iHnghes et al.l 1 1 9981: iBarger et al.l Il99^ were car¬ 
ried out with the Sn bmillimeter Gammo n User Bolome¬ 
ter Array (SCUBA: iHolland et al.l [19^ at 850/im, the 
existence of a population of high-redshift dusty star¬ 
forming galaxies (DSEGs) has been established, expand¬ 
ing our understanding of galaxy formation and evolution. 
Thes e massive galaxies (M:^. iHainline et al.l 

I 2 QIII : iMichalowski et aP I 2 Q 12 I: iSimpson et al.l l2QI4f ) are 
extremely bright at snbmillimeter wavelengths but faint 
in the optical due to extinction and hence are also 
known as “snbmillimeter galaxies” (SMGs). Their spec¬ 
tral energy distributions (SEDs) are dominated by rest- 
frame far-infrared luminosities in excess o f IO^^Lq (e.g., 
iCoppin et al.l I 2 QQ 8 I: iMagnelli et al.l I 2 QI 2 D . These dusty 
luminous galaxies are producing stars at typically in- 
tense star formation rates of > 100 — 1 000 M 0 yr“^ (e.g., 
iNeri et aLll2QQ3l : IChaprnan et al.ll2QQ5f ). The high SERs 
can be explained by either large reservoirs of molecu¬ 
lar gas or a bo osted star formation efficiency induced by 
major mergers (iNaravanan et al.ll2QQ9l: ll 


lEn et al.ll2Q13l : IB ether min et al. l2Q15bl : 


Ingel et a 


Dve et al 


Mia 

mm . 


The dominant mechanism is still a topic of debate. 

Previous studies of their stellar masses and prodi¬ 
gious star formation rates have been limited to the 
most luminous systems due to the challenge of obtain¬ 
ing redshifts and detailed follow-up observations of these 
distant, dns t-obscnred objects. T he South Pole Tele¬ 
scope (SPT: [Carlstrom et al.ll2Qll| ). which surveyed 2500 
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deg^ at 1.4, 2, and 3mm with high spatial resolution 
(FWHM ^1.0'), has helped us overcome these difficul¬ 
ties by discovering ^100 high redshift, strongly gravita¬ 
tiona lly lensed DSFGs (jVieira et al.ll2QlQl: iMocann et al.l 
I2Q13F The leasing origin of the SPT DSFGs was con¬ 
firmed by the A tacama Large Millimeter/submillimete r 
Array (ALMA) (|Hezaveh et al.|[2QTl:IVieira et al.ll2Q13[) . 
The ALMA submillimeter imaging is used to model the 
l ensing geometry yiel ding magnification factors up to 30 
(jHezaveh et al.lI2QT^ Spilker et al. in prep). Strong 
gravitational lensing enables detailed analyses of these 
DSFGs but also complicates recovery of the properties of 
the sources because emission from the background DSFG 
is blended and confused with that from the foreground 
lensing galaxy. In the FIR/submillimeter regime, the 
background emission dominates, but in the rest-frame 
optical/near-IR where we probe the established stellar 
population, we must de-blend the background DSFG 
from the foreground lens. 

The primary goal of this paper is to measure stellar 
masses and hence determine the specific SFRs (sSFRs), 
which are a critical diagnostic for understanding the na¬ 
ture of these galaxies. We develop and implement an im¬ 
age de-blending technique with which we can recover the 
stellar mass of the lensed DSFG. With the lens model, 
we are able to derive the intrinsic properties of the SPT 
DSFGs (e.g., stellar mass, luminosity, and star formation 
rate) and place these galaxies in the context of the evo¬ 
lution of the high-redshift galaxy population as a whole. 
We present six SPT DSFGs all of which have been sys¬ 
tematically followed up with HST, Spitzer^ APEX, and 
Herschel and have robust lens models derived from high 
resolution ALMA imaging. 

This paper is organized as follows. In Section [2l we 
present our multi-wavelength observations. Section [3] de¬ 
scribes the de-blending technique and results for each 
source. In Section IH we describe the lens model. In Sec¬ 
tion [5l we present the stellar masses and star formation 
rates derived from SED fitting. We discuss the SER - 
relation and sSERs in Section [6l A summary of the fea¬ 
sibility analysis of the de-blending technique and conclu¬ 
sions are presented in Section 0 Throughout this work, 
we use a standard cosmological model with Hq = 70 km 
s“^Mpc“^, Qm = 0.3, and Qa = 0.7. 

2. OBSERVATIONS 

The six SPT DSEGs presented in this paper (Tables [U 
and [2]) are drawn from a parent sample of 50 SPT DS¬ 
EGs with Spitzer imaging, and constitute the best cases 
among the subset (13 sources ) with HST data and lens 
models (Sec. |4]) derived from ALMA imaging. Eig. [T] 
shows the HST^ Spitzer^ and ALMA data for the six 
sources. 

2.1. Spitzer/IRAC 

A total of 50 SPT DSEGs were followed up by two 
Spitzer Space Telescope programs (PID 60194 and PID 
80221; PI Vieira) and a joint Spitzer/HST program 
(PID 1 0094; PI Vieira) u sing the Infrared Array Camera 
(IRAC: [Fazio et al1l2004l ). For the Spitzer programs, im¬ 
ages were taken at 3.6/im and 4.5/im with 36 or 32 (dither 
pattern) x 100 second (frame time) exposures and 12 
(dither pattern) x30 (frame time) exposures, respec¬ 
tively. For the joint program, we took the 100s frames 


in both bands with 36/72/108 dither positions. We sum¬ 
marize the data for each source in Table [H The basic cal¬ 
ibrated data, pre-processed using the standard pipeline 
by the Spitzer Science Center, were resampled and com¬ 
bined into a mosaic image utilizing the MOPEX soft- 
ware package (iMakoyoz fc Marleaull2005[ ) and IRACproc 
([Schuster et al.ll2006f ). The IRAC mosaics have a resam¬ 
pled pixel scale of 0.6"/pixel and an angular resolution 
of ^ 1.7". The IRAC data are the primary probes of 
the stellar component of the background DSFGs. How¬ 
ever, due to the close proximity of the foreground lens¬ 
ing galaxy, emission from the lensed DSFG is mixed with 
that of the foreground lens and needs to be de-blended. 
The de-blending technique, which will be described in 
SecHl is a critical component of this project, making 
extracting fiuxes from extremely distant sub-millimeter 
sources possible. 

2.2. HST/WFC3-IR 

A sample of 18 SPT sources were observed with 
the Wide Field Camera 3 (WFC3) on board the Hub¬ 
ble Space Telescope under program 12659 (PI: Vieira). 
Each target was assigned one complete orbit split be¬ 
tween El low and F160W filters in the infrared chan¬ 
nel. Another six sources were followed up in the joint 
Spitzer/HST (PID 13614; PL Vieira) at F140W 

with one orbit assigned per source. Dithering was imple¬ 
mented for maximum resolution. The data were reduced 
using the standard HST pipeline. The pixel size of the 
WFC3 images is 0.128". Detailed information on the 
WFC3 data pro ducts can be found in t he WFC3 instru¬ 
ment handbook (jPressel fc et al.ll2015[ ). 

2.3. ALMA imaging 

The gravitationally lensed origin for the SPT sources 
was confirmed by ALMA under a Cycle 0 program 
(2011.0.00958.S; PI: D. Marrone). We imaged 47 SPT 
sources with ALMA at 870/im with the dual-polarization 
Band 7 (275-373 GHz) receivers in both compact and 
extended array configurations , reaching ^ O . b" res oln- 
tion. We refer the reader to iHezaveh et al ] (pnTl for 
details of the ALMA observations. The visibilities from 
each configuration were concatenated and imaged using 
uniform weighting to emphasize resolved structure. The 
combined ALMA images are used in this project to indi¬ 
cate position and structure of the DSFGs. 

2.4. FIR data 

2.4.1. APEX/LABOCA 

The SPT sources were imaged at 870//m with the Large 
APEX Bolometer CAmera (LABOCA) at APEX under 
the program M-085.F-0008 (2010, PI:Weiss). LABOCA 
is a 295-element bolometer array (Siringo et al. 2009 ) 
with an 11.4' field-of-view and an angular resolution of 
19.7" (FWHM). The central frequency of LABOCA is 
345GHz (870/im) with a passband FWHM of ^ 60GHz. 
Observations were carried out under good weather con¬ 
ditions (typical PWV of 0.9mm ranging between 0.3 and 
1.5mm). The data reduction wa s performed in the same 
manner as in iGreve et a, 1.1 (1201 2h . 

2.4.2. Herschel/SPIRE 
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Figure 1. 20'' x 20" cutouts of the six SPT DSFGs showing the HST/WFC3 (grey), Spitzer/IRAC (blue contours), and ALMA Band 7 
(red contours) data. 


Table 1 

Summary of Spitzer/IRAC and HST/WFC3 data 


SPT Source Short Name 


SPT-S J010312-4538.8 
SPT-S J034640-5204.9 
SPT-S J053816-5030.8 
SPT-S J213404-5013.2 
SPT-S J214654-5507.8 
SPT-S J214720-5035.9 


SPT0103-45 

SPT0346-52 

SPT0538-50 

SPT2134-50 

SPT2146-55 

SPT2147-50 


Filters (exposure time) 

FllOW (1312s), F160W (1412s), 3.6/im (3600s), (7200s) 

FllOW (1312s), F160W (1412s), 3.6/rm (3600s), 4.5/rm (10800s) 
FllOW (1312s), F160W (1412s), 3.6/rm (3600s), 4.5/rm (360s) 
FllOW (1312s), F160W (1412s), 3.6fiin (3200s), Abfim (360s) 
F140W (2812s), 3.6p.m (3600s), (3600s) 

F140W (2812s), 3.6/im (3600s), 4.5/rm (3600s) 


Spitzer PIP HST PIP 


10094 

10094 

60194 

80221 

10094 

10094 


12659 

12659 

12659 

12659 

13614 

13614 


The SPIRE (| Griffin et alJ I2Q1Q[ ) data were obtained 
in the programs OTl_jvieiraA and OT2_jvieiraA>^ and 
were observed simultaneously at 250/im, 350/im and 
500/im. The maps were produced using the standard 
reduction pipeline HIPE v9.0. The fluxes were extracted 
by fitting a Gaussian to the source and reading out the 
maximum value. The fluxes have been corrected for pix- 
elization as described in t he SPIRE Observers’ Manual 
(|Valtchanov fc et al]l2Q14D . 

2.4.3. Herschel/PACS 

The PAGS (|Griffin et al.l |2Q1Q[) data at lOO/im and 
160/im were obtained in the programs OTl_jvieiraA^ 
OTlAmarroneA and OT2_jvieiraJ). Each PAGS map 
has been co-added, weighted by coverage. We performed 
aperture photometry with the aperture sizes fixed to 1" 
for the 100/im map and 10" for the 160/im map. Details 
of the EIR data observations, reduction, and photometry 
will be given in Strandet et al. (in preparation). 

2.5. ALMA redshift survey 

Using ALMA, we conducted a blind redshift survey in 
the 3mm atmospheric transmission window for a subset 
of the SPT sources. The lines detected are identified as 
redshifted emission lines of ^^GO, ^^GO, GI, H 2 O, and 
H 2 O+. The details of the spectral line observations and 
the redshift distribution of dusty sty- for ming galaxies 
by SPT can be found in IWeifi et al.l (j2013f ). We present 


in Table [2] the redshifts for the sources analyzed in this 
paper. 

3. DE-BLENDING TEGHNIQUE 
3.1. Modeling the lenses in HST 

To construct models of the surface brightness profiles 
for each le nsing galaxy, we perform profile fitting using 
GALEIT (jPeng et al.l DOlOD . which can simultaneously 
fit profiles to multi ple galaxies . We start with a 
single Sersic profile (jSersid Il963f ) and then add other 
components to generate better models as needed. We 
create the models using EllOW filter images when 
available, as EllOW contains less emission from the 
background DSEG than E160W and therefore enables 
better modeling of the foreground lens. 

Here we present the modeling results for the six 
systems in our DSEG sample. S PTQ5 38-5Q, which has 
been modeled in iBothwell et al.l (j2Q13f ). serves as a test 
case for our analysis. Eor the non-detections, we place 
upper limits on the flux densities to constrain the SED 
of the background source. To derive photometric upper 
limits, we performed random aperture photometry on 
the background using a radius for each source that 
encloses the 3cr ALMA contours. The detections and 3cr 
upper limits are presented in Table [3l 

3.1.1. SPT0538-50 
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Table 2 

Lens and source parameters 


SPT source name 

ZL 

zs 

n 

Sersic 

(n, re, mag) 

ExpDisk {n = 1) 
(re, mag) 

PSF 

(mag) 

SPT0103-45^ 

0.740 

3.092 

5.34 ± 0.11 

(1.9, 0.3", 21.84) 

— 

23.89 

SPT0346-5:P 

-1.1 

5.656 

5.57 ± 0.12 

(2.3, 0.6", 24.58) 

— 

— 

SPT0538-5(Jal 

0.404 

2.786 

20.12 ± 1.81 

(3.2, 0.7", 18.58) 

(2.4", 19.74) 

20.97 

SPT2134-5(P 

0.776 

2.780 

21.00 ± 2.42 

4 . 9 , 0.6", 20.17) 

— 

21.88 

SPT2146-55l> 

— 

4.567 

6.65 ± 0.41 

(3.1, 0.1", 23.77) 

(0.6", 22.09) 

— 

SPT2147-5(P 

0.845 

3.761 

6.55 ± 0.42 

(2.7, 0.4", 20.02) 

5 . 3 ", 20.80) 

— 


Note. — We adopt the source redshift for SPT0538-50 obtained with ATCA CO(l-O) spectroscopy bv lAravena et al.l (I2013I L We use 
the lensing magnification factors in Spilker et al. (in preparation). The lens parameters (e.g., Sersic index, effective radius) derived from 
FllOW or F140W images using GALFIT are held fixed during the de-blending process. We list the parameters for each component: the 
Sersic profile, the ExpDisk (i.e., exponential disk), and the PSF point source component. 

^Magnitudes (AB) measured in the FllOW filter. 

^Magnitudes (AB) measured in the F140W filter. 


We first simultaneously fit the central lensing galaxy 
and 9 neighboring objects with a single Sersic profile for 
each of them as a first pass. The resultant best-fit model 
for the central lens has a Sersic index (i.e., the shape 
parameter controlling the central concentration of the 
profile) of n = 8.8 and large residuals with significant 
positive and negative counts left over in the center. We 
improve the fit by using a three-component model for 
the central lensing galaxy, consisting of a Sersic profile, 
plus an exponential disk, plus a point source (SEP). A 
Sersic bulge and an exponential disk are insufficient and 
the point source is needed to minimize the residuals. The 
SEP model yields reasonable residuals with n = 3.2. We 
present the fitting result of the SEP model in Eig. [2l 
The top panels of Eig. [2] show the original WEC3 image 
(top left), the best-fit model by GALEIT (top middle), 
and the residuals leftover after subtracting the best-fit 
model from the original image (top right). We see no 
excess emission associated with the structure indicated 
by the ALMA contours, suggesting the fluxes extracted 
from the WEC3 EllOW image are only from the fore¬ 
ground lensing galaxy. Adding additional components 
does not substantially improve the fit, and is not physi¬ 
cally motivated. We therefore continue the de-blending 
pr ocess with the three-c omponent SEP model. Note that 
in iBothwell et al.l (j2Q13f ) the lensing galaxy was modeled 
using a co-add between the EllOW and E160W bands 
while we use the bluest band throughout the paper. 

3.1.2. SPT0103-45 

The source is best fit with a Sersic model {n = 1.9) 
and a point source. The obvious “arc”-like structure 
(Eig. [3] top left) is not the lensing arc associated with 
the background DSEG, which is indicated by the ALMA 
contours. We treat that as two elongated Sersic-profile 
galaxies that are subtracted off from the original image. 
The top middle panel shows the best-fit model by GAL¬ 
EIT and the overall residuals after the subtraction are 
shown in the top right panel. 

3.1.3. SPT0346-52 

SPT0346-52 has the highest redshift {z = 5.7) among 
all the sources in this work. We fit four galaxies in the 
field with a single Sersic profile {n = 2.3 for the lensing 
galaxy) by GALEIT, shown in EigIl](top panels), yield¬ 
ing a clean residual. Given its high redshift, it is not 
surprising that SPT0346-52 is undetected in the EllOW 


band. We note that the arcs are well-separated from 
the central lens with a diameter of 2.2", allowing a de¬ 
coupling of the background and foreground emission in 
IRAG. 

3.1.4. SPT2134-50 

We simultaneously fit the lensing galaxy and 4 neigh¬ 
bors in the field (EigUltop panels). A point source was 
added to the central lens in addition to the Sersic profile 
(n = 1.9). We see no trace of the emission suggested by 
the sub-millimeter ALMA contours in the residual image 
(top right). Note that the lensing structure indicated by 
the ALMA contours is marginally resolved with a separa¬ 
tion of 1.3" in diameter, smaller than the separations in 
SPT0538-50 of 1.8" at a similar redshift and SPT0103-45 
of 2" at even higher redshift. 

3.1.5. SPT2146-55 

SPT2146-55 can be fit with a Sersic {n = 3.1) plus an 
exponential disk in the E140W image. There is excess 
emission left over in the residual image, matching the 
ALMA contours. We find a detection of the background 
DSEG in this system at E140W with a flux density of 
(1.5 ± 0.4) X 10“^ mJy. A simple Sersic model also re¬ 
veals the background source, albeit with worse central 
residuals. 

3.1.6. SPT2147-50 

The lensing galaxy of SPT2147-50 is fit with a three- 
component SEP model of which the Sersic index is n 
= 2.7. There is no residual associated with the ALMA 
contours. 


3.2. De-hlending the source/lens in IRAC 

With the modeling of the lenses in the high-resolution 
HST images, we are able to disentangle the DSEG 
emission from the foreground emission in IRAG bands. 
The first attempt to disentangle the contributions from 
the foreground lens and the background source in low- 
resolution Spitzer/IRAC images of strongly lensed DS- 
EGs was made by Hopwood et al. (2011). Similar 
de-blending techniques with additional r efinements were 
implemented in sulpseqne nt works (e.g., iBothwell et al.l 
[201^ iNegrello et al.l[20Tl) . 

We utilize the code Python Galaxy Eitter (PyG- 
Eit; iMancone et al.|[2QT^ which is designed to measure 
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Figure 2. 23" X 23" HST diiad IRAC cutouts of SPT0538-50 showing the de-blending technique. Each image is oriented such that North 
is up and East is to the left. In each panel, the green contours show ALMA sub-millimeter continuum emission, indicating the position 
and structure of the background dusty star-forming galaxy. Contours are in steps of Scr, starting at 3a. The top row shows modeling of the 
lens in HST WFC3 FllOW. Top Left: The original HST/WFC3 image of SPT0538-50 in FllOW filter. Top Middle: The best-fit GALFIT 
model of SPT0538-50. Top Right: The residual image after subtracting the best-fit model from the original HST image. The residual flux 
of the knot south-east away from the center is 0.3%. 

The middle row and bottom row show de-blending of the source/lens in IRAC at 3.6/i,m and d.S/xm, respectively. Left: The original 
Spitzer/IRAC 3.6/i,m/4.5)nm image of SPT0538-50. Middle: The GALFIT model convolved with the IRAC 3.6/i,m/4.5/rm PSF. Right: The 
residual image. The Einstein ring is recovered in IRAC after the foreground lens is subtracted off. 


PSF (point-spread function)-matched photometry from 
images with disparate pixel scales and PSF sizes. It 
takes models generated from a higher resolution image, 
HST in this case, and fits blended sources in crowded, 
low resolution images (i.e., Spitzer/IRAC). It should be 
noted that we make the assumption that the shape of 
the light profile of foreground galaxies remains the same 
from WFC3 1.1/im to IRAC 3.6/im, and 4.5/im, i.e. mor¬ 
phological /^-corrections are negligible. Before fitting the 
IRAC images by performing a minimization, the HST 
models are convolved with the IRAC PSF, which is em¬ 
pirically extracted by stacking at least five bright and 
isolated stars in the field. The quality of the IRAC PSF 
is critical to de-blending to avoid artifacts introduced by 


the PSF. We test the quality of the IRAC PSF by ran¬ 
domly selecting stars in the same field and fitting them 
with the same PSF used in de-blending. We examine the 
star-subtracted residuals to assure the PSF is acceptable. 
A few examples of such fit and subtraction are shown in 
Fig. [HI The residual flux from these stars is less than 1% 
of the original flux in all cases. 

We optimize the fitting by masking out the region 
where the DSFG emission is expected based upon the 
ALMA data, thereby fitting only the foreground lens. 
During the minimization only the positions and fluxes 
of the objects are left as free parameters. All other Sersic 
parameters (effective radius, Sersic index, aspect ratio, 
and position angle) are held fixed. The positions are re- 













6 


J. Ma, et al. 



Figure 3. 6 " x 6" HST and IRAC cutouts of SPT0103-45 showing the de-blending technique. Subfigures are the same as in Figure O 


stricted to small shifts (typically less than a pixel). We 
examine the residuals to see if there is excess emission 
associated with the ALMA contours. 

3 . 2 . 1 . SPT0538-50 

The middle row and bottom row in Fig. [2] show the de- 
convolution process just described at 3.6/im and 4.5//m 
adopting the SEP model. The best-fit IRAC model 
{Middle)^ constructed by convolving the iJ^T/WFCS 
model with the IRAC PSF, is removed from the orig¬ 
inal 3.6/im/4.5//m image {Left) to produce the residual 
image on the right. The Einstein ring is clearly recovered 
in both IRAC bands, closely matching the ALMA con¬ 
tours in green. We perform aperture photometry on the 
IRAC residual images, using an annular aperture large 
enough to enclose the 3cr ALMA contour without in¬ 
cluding too much sky flux. The resultant sky-subtracted 
flux densities (Table [3]) at 3.6/im and 4.5/im are Ss,q = 
26.2 ± 6.7/iJy and 5'4.5 = 53.7 ± 7.0/iJy. We note that 
these are consistent with the values of Ss,e = 22 ± 5/iJy 


and 5'4.5 = 47 ± 8/iJy derived bv iBothwell et al] (j2Q13f ). 
The uncertainties include the photometric uncertainties 
as well as the uncertainty associated with the de-blending 
process where the residual fluxes can vary due to a dif¬ 
ferent set of PyGFit input parameters. We find that 
the maximum allowed positional shift during the fit¬ 
ting dominates the uncertainties in the resultant residual 
structures, and therefore the residual fluxes particularly 
for simultaneous fitting of the source consisting of multi¬ 
components. The above best-fit model is determined by 
minimization. We apply the Ay^ technique to de¬ 
termine the Icr error by varying the input positions in 
PyGFit. We change the allowed maximum positional 
shift during fitting from 0 to 0.6" (1 pixel) until y^ in¬ 
creases by 2.3. 

3 . 2 . 2 . SPT0103-45 

We apply the analogous de-blending technique to 
SPT0103-45 and it is detected in both bands, with excess 
emission coincident with the ALMA contours. We put 
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10" X 10" HST and IRAC cutouts of SPT0346-52 showing the de-blending technique. Subfigures are the same as in Figure O 


Figure 4. 

an annular aperture on the residual images, extracting 
the flux densities at two bands of 83 , 6 = 5.1 ± l.SyuJy and 
5'4.5= 12.8 ± 1.4/iJy. The errors are determined in the 
similar way. 

3 . 2 . 3 . SPT0346-52 

The middle row and bottom row in Fig. |4] show the 
results of the de-convolution process for SPT0346-52 in 
IRAC bands. The 3.6/im and 4.5/im images are sub¬ 
tracted clean. SPT0346-52 is a non-detection in IRAC 
with 3 hour exposure at 4.5 /im even with the help of 
strong gravitational lensing. Formally we measure the 
flux to be 0.5 ± 0.8 /iJy at 3.6/im and 0.3 ± 1.2 /iJy at 
4.5/im. We use the non-detections in both IRAC bands 
to put upper limits on the fluxes from SPT0346-52 based 
on the sky flux distribution (random aperture) in the 
residual maps. We place a 3.6/im band 3cr upper limit of 
2.4 /iJy and a 4.5/im band 3cr upper limit of 3.6 /iJy . 

3 . 2 . 4 . SPT2134-50 


Examining the IRAC residual maps (Fig. [5]) of 
SPT2134-50, we find the excess emission coincident with 
the ALMA contours. However, due to the small separa¬ 
tion of the lensing arcs, imperfections in the HST model 
could also contribute to the residual in IRAC given the 
relatively large IRAC PSF. In this case, we failed to ex¬ 
tract robust fluxes from the residual because the separa¬ 
tion was not large enough (^ 1.7" in diameter) to avoid 
confusion. We instead extracted the fluxes as 3cr upper 
limits of < 20.0 /iJy at 3.6/im and < 32.9 /iJy at 4.5/im, 
providing a loose constraint on the stellar mass. 

3 . 2 . 5 . SPT2146-55 

With the detection in HST F140W, we expected to 
detect SPT2146-55 in the IRAC bands. However, the 1 
hour IRAC data do not provide clear detections, given 
that the residuals do not match the ALMA contours and 
are not clearly related to the background source. The 
measured fluxes are 0.5 ±1.4 /iJy at 3.6/im and 0.4 ± 
2.3 /iJy at 4.5/im. We instead adopt 3cr upper limits 
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9" X 9" HST and IRAC cutouts of SPT2134-50 showing the de-blending technique. Subfigures are the same as in Figure [2] 


Figure 5. 

estimated from the background sky using a ring aperture: 
< 5.2 yuJy at 3.6 jam and < 6.9 //Jy at 4.5 jam. 

3.2.6. SPT2147-50 

After the de-blending process for SPT2147-50, there 
are three blobs coincident with the ALMA contours. 
We use separate apertures corresponding to the separate 
ALMA contours to extract fluxes instead of an annular 
aperture to avoid the over-subtracted regions. We obtain 
flux densities of 2.0 ± 0.6 /iJy at 3.6 jam and 3.5 ± 0.4 
/iJy at 4.5 /im. 

4. LENS MODELS 

The lens modeling technique, based o n our Cycle 0 
ALMA Band 7 imaging, is described by iHezaveh et al.l 
(j2013f ). In this paper, six gravitationally lensed SPT 
sources, including SPT0538-50 and SPT0346-52, are dis¬ 
cussed in detail. To summarize, the modeling is per¬ 
formed in the visibility plane to properly compare with 
data as the ALMA interferometer samples the Fourier 


transform of the sky brightness distribution across a 2-D 
spatial frequencies instead of directly imaging the emis¬ 
sion. The source is modeled as a symmetric Gaussian or 
a Sersic light profile while the lens is assumed to be a Sin¬ 
gular Isothermal Ellipsoid (SIE). Minimizing between 
data and model visibilities determines the best-fit. The 
lensing magnification factor, /i, is defined as the ratio of 
the total lensed to unlensed flux. The lens models allow 
us to derive intrinsic properties of the lensed galaxies by 
simply dividing the observed properties by fi. We present 
the resultant magnifications for the sources discussed in 
this work in Table [2l We adopt the latest magnification 
factors for all the six sources derived in the same fash¬ 
ion in Spilker et al. (in preparation). When deriving 
stellar masses and SERs, we correct for the lensing mag¬ 
nifications and propagate the corresponding errors. We 
con sider the issue of differential magnification in Section 

[Ql 

5. SED FITTING 
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Figure 6. 


10" X 10" HST and IRAC cutouts of SPT2146-55 showing the de-blending technique. Subfigures are the same as in Figure O 


Table 3 

Multi-band observed fiux densities in mJy 


Observed Wavelength A 

SPT0103-45 

SPT0346-52 

SPT0538-50 

SPT2134-50 

SPT2146-55 

SPT2147-50 

HST/WFC3 1.1/rm 
HST/WFC3 1.4/im 

< 5.8 X 10-4 

< 3.8 X 10-4 

< 0.0022 

< 0.0012 

(1.5 ±0.4) X 10-3 

< 7.3 X 10-4 

HST/WFC3 1.6/im 

< 0.0026 

< 9.1 X 10-4 

< 0.0033 

< 0.0016 



Spitzer/IRAC 3.6/i,m 

0.0051 ± 0.0013 

< 0.0024 

0.0262 ± 0.0053 

< 0.0200 

< 0.0052 

0.0020 ± 0.0006 

Spitzer/IRAC A.bfim 

0.0128 ± 0.0014 

< 0.0036 

0.0537 ± 0.0070 

< 0.0329 

< 0.0069 

0.0035 ± 0.0004 

Herschel/PACSlOOiim 

< 13 

< 6 

31 ±2 

49 ±3 

< 8 

9±2 

Herschel/PACSlGOfim. 

< 47 

33 ±9 

141 ± 15 

196 ± 22 

< 29 

< 28 

Herschel/SPIRF 250/rm 

133 ± 14 

122 ±11 

326 ± 23 

350 ± 25 

65± 13 

72 ±9 

Herschel/SPIRF 350)am 

213 ±16 

181 ± 14 

396 ± 38 

332 ± 23 

69± 13 

115 ±10 

Herschel/SPIRF 500/xm 

232 ± 17 

204 ± 15 

325 ± 24 

269 ± 19 

83± 10 

121 ± 11 

APEX/LABOCA STO/xm 

125 ±6 

131 ±8 

125 ±7 

101 ±7 

55 ±4 

61 ±5 

SPT 1.4mm 

39.1 ±7.0 

46.0 ±6.8 

28 ±4.6 

22.8 ±4.6 

17.8 ±3.9 

20.3 ±4.6 

SPT 2.0mm 

8.8 ± 1.4 

11.6 ±1.3 

8.5 ± 1.4 

6.0 ± 1.3 

5.2 ± 1.3 

5.8 ± 1.3 


Note. — All the HST and IRAC photometry is derived in this work. The quoted uncertainties include photometric uncertainties plus 
uncertainties due to the de-blending process. For the non-detections, the fiux density upper limits are given at 3a. The FIR photometry 
is from Strandet et al. in prep except SPT0538-50 for which we use the photometry bv IBothwell et aTI 120131 b 
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Figure 7. 1?>" x 1?>" HST and IRAC cutouts of SPT2147-50 showing the de-blending technique. Subfigures are the same as in Figure O 
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Figure 8. Examples of stars in the field of SPT0538-50 at 4.5 p,m fit with the same IRAC PSF used in de-blending. Top: Original images. 
Bottom: Residual images after the original images subtracted by best-fit PSF. 


To derive properties of these dusty star-forming galax¬ 
ies such as stellar mass, infrared luminosity {Lir from 
8-lQQQ /im), and SFR , we utilize the c ode CIGALE 
(jBurgardla et al.l l2QQ5l : iNoll et aH l2QQ9[Fn to fit our 
photometric data from near-IR to far-IR. For the fit¬ 
ting, the code generates FUV to FIR SEDs consisting 
of dust-attenuated complex stellar population models, 
IR dust emission models, and spectral line templates. 
It allows self-consistent treatment of the stellar com¬ 
ponent, sampled by our NIR data and the re-emitted 
dus t compone nt, sampled by the FIR data. We adopt 
the iMarastonI f2 Q Q5; M O5) stellar population synthesis 
model, a lKrouPc] (|2QQI[) initial mass function (IMF), and 
a solar metallicity as our baseline. 


5.1. Star formation histories 

One- or two-component star formation history (SFH) 
models with exponentially increasing or decreasing SFR 
or continuous constant SFR are possible to define the 
star formation history in CIGALE. Commonly assumed 
star formation histories in the literature are: (I) a sin¬ 
gle stellar population w ith an exponenti a lly d e creas¬ 
ing SFR (r model ; e.g., Illbert et al.l l2QQ9l I2QIQL I2QI3I: 
iMuzzin et al.ll2QI3r). (2) an exponentially increas ing SFR 
(e.g.. lMaraston et^l2QIQl:lPapovich et al.ll2QIIf ). (3) de¬ 
layed r model (e.g. Ida Cunha et al.ll2QI5f ). (4) an instan¬ 
taneous b urst of star formation or constant st ar forma¬ 
tion (e.g. lKennicut0lI998l: iMurphv et al.ll2QIIf ). and (5) 
a model with two stellar populations: a young stellar 
population with a constant SFR on top of an old stellar 
population with an exponentially d eclining SFR (e.g., 
iPapovich et al .1120011: iLee et al~]l2QQ9[) . 


5.1.1. Single-component SFH model 

The most commonly used SFH model is the exponen¬ 
tial or r-model with the simple form of SFR oc 
This model naturally arises in scenarios where the SFR 
is directly proportional to the gas density in a closed- 
box model (jSchmidtllI959[ ). There are two parameters, 

http://cigale.lam.fr 


the e-folding time r and age which we consider as free pa¬ 
rameters. If r >> age, the r-model transforms to a single 
burst star formation history with a constant rate. This 
simple model is widely adopted in the lite rature, facili¬ 
tating comparison with previous work (e.g..lKarim et al.l 
I2QIII: Heinis et al.ll2Ql4 iBethermin et al.ll2QI5b[) . 

5.1.2. Two-component SFH model 

We also investigate the two-component SFH model 
because real systems are expecte d to experience mul¬ 
tiple epis o des o f star formation (jPapovich et al.l I2QQII: 
iLee et all I2QQ9I ). Moreover, the double SFHs are 
l ikely to better repro duce the true stellar mass 
(jMichalowski et aP l2QI4l ) . The two stellar population 
components are connected by the burst strength/fraction 
fysp (i.e., the mass fraction of young stellar population in 
the model). We first treat the ages of both the young and 
old stellar populations and the burst fraction as well as 
Av as free parameters. The derived instantaneous SFRs 
and stellar masses are poorly constrained. These two 
parameters, the age of the young stellar population and 
the burst fraction, are the dominant ones that drive the 
resultant stellar mass. The dependence of derived stel¬ 
lar mass and SFR on varying ages and burst fractions 
is investigated. We show the case of SPT0538-50 as a 
demonstration in Figure [9l The other SPT DSFGs ex¬ 
hibit similar trends. We allow the age of the young stellar 
population to vary from 10 to 100 Myr in steps of 10 Myr 
(commonly adopted starburst lengths ^10-100 Myr, e.g., 
Michalowski et alll20I0l: IMurphv et al.l [ 2 MTI: iKennicu^ 
I998f ) . We permit the burst fraction to vary from 0 to I, 
noting that burst fractions of < O.I are ruled out at high 
confidence. We generate a grid of models of different in¬ 
put fysp and tysp. The resulting SEDs are equally well 
fit in terms of matching the IRAC and FIR data points 
and being consistent with the HST upper limits. As a 
result, the stellar mass increases by ^0.8 dex with an 
increasing age from 10 to 100 Myr. It is also subject to 
the burst fraction at a fixed age: drops by ^ 0.8 dex 

from fysp = O.I to I. The instantaneous SFR is much 
less dependent on fysp than it is on tysp. This suggests 
that we do not have enough leverage to constrain these 
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Figure 9. The dependence of derived stellar mass on the age of 
the young stellar population and the burst fraction for SPT0538- 
50. 

key parameters and the stellar mass heavily depends on 
the assumptions made for the two parameters. 

We demonstrate in the SFR - M* plane (Fig. [TQ|) that 
SFRs modeled by the single-component SFH are better 
constrained than those modeled by the double-SFH. We 
adopt the simplest single-component SFH model here¬ 
after since we do not have enough data to constrain a 
more complex SFH. We present the SED fitting results 
in the next section with the r-model adopted. 

5.2. SED fitting results: Stellar mass and SFR 

We run CIGALE, with the input parameters listed in 
Table IH to find the best-fit SED and inferred stellar 
mass, IR luminosity, and instantaneous star formation 
rate. Expectation values and standard deviations for the 
output galaxy properties are derived from the probability 
distribution function of the parameter-bin-specific best- 
fit models. Eig. [TT] shows the individual best-fit SEDs 
for the six SPT DSEGs in the observed frame. We also 
obtain an intrinsic (i.e., magnification corrected) aver¬ 
age SPT SED (in black in Eig. [12]) by taking the median 
value of the six best-fit SEDs at each rest-frame wave¬ 
length. We report stellar masses, infrared luminosities, 
SERs and associated uncertainties in Table El The EIR 
part of the SED is well constrained, thus the IR lumi¬ 
nosity has little dependence on varying input parame¬ 
ters. The optical/UV part is less constrained, therefore 
the derived stellar mass depends strongly on the input 
parameters. The dust extinction is a free parameter in 
the fitting and we find a significant attenuation by dust 
{Ay ^ 4 — 5) is required to match the IRAC photometry 


and the EIR data points, which i s consistent with other 
SMGs (e.g.. iNegrello et ahlDOlOl: iHopwood et al.l I 2 Q 11 I: 

iHainline et al.ll^lll: iMichalowski et al.ll2QlQf ). The de¬ 

rived ages are all much smaller than the corresponding 
e-folding times suggesting that the star formation history 
essentially can be described by a single extended burst 
star formation. The best-fit burst ages are 20 - 80 Myr. 
However, we note that there can exist older stellar pop- 
nlations that are hidden by the bright y oung population 
(jMaraston et al.ll 2010 l: iBnat et ^l2014| ). 

Our stellar mass estimates for the six SPT DSEGs span 
more than an order of magnitude, with a median value 
around 5 x IO^^Mq, for those detected. The median 
stellar mass is lower t han the median values of a well- 
studied SMG sam ple (jChapman et al.l l2005f ) found by 
IMichalowski et ahl ( 2010 ; 3.5 x IO ^Mq) and the recent 
ALESS sample bv ida Cnnha et al] (2015; 8.9 x IO^^Mq). 
These two studies perform multi-band SED fitting from 
rest-frame UV to EIR as we do. However, these two sam¬ 
ples are biased towards higher stellar masses while ours 
are biased towards leasing configuration. 

All of the SPT DSEGs present prodigious star forma¬ 
tion at rates > 500 M 0 yr“^. SPT0346-52 stands out as 
having a SER of ^4800 Mqyt~^ only ^ 1 billion years af¬ 
ter the Big Bang, among the highest SERs at any epoch. 
Given that the Lm is the best measured quantity and 
is widely used as an SER indicator, we also calculate 
the SE Rs u sing Ljr conversi on recipes of iMnrphv et al.l 
(j2011f ) and iKennicnttl (jl998f ) with the common IME of 
Kroupa and continuous bursts of age 10-100 Myr. When 
converting SER from Salpeter IME to Kroupa IME, we 
divide the conversion coefficient by 1.53. Both estimates 
are presented in Table El All of the SERs derived by 
GIGALE SED fitting are consistent with those by the 
conversion recipes. 

5.2.1. Systematics 

A scatter of ^ 0.3 dex around the correct value is 
an intrinsic uncertainty in determining stellar masses 
from broad-band phot ometry through SED fitting 
(jMichalowski et aT]l2014f ). Eor our SPT DSEGs, we have 
checked that stellar masses derived using the Bruzual & 
Chariot (2003; BC03) models are at most 0.1 dex higher 
than masses derived using the M05 mo del A two- 

component SEH on average yields stellar masses ^0.15 
dex higher than a single-component SEH. These system¬ 
atic uncertainties are subdominant compared to the sta¬ 
tistical uncertainties (^ 0.4 dex) for our data. Our un¬ 
certainties are dominated by the limited number of pho¬ 
tometric data points at optical/near-IR wavelengths. 

6 . DISCUSSION: SFR-M* RELATION — A COMPARISON 
WITH STAR-FORMING MAIN SEQUENCE 

6 . 1 . SFR-M^ relation 

Galaxy surveys at low and high redshifts have shown 
that star-forming galaxies form a power-law relation be¬ 
tween their star formation rate and stellar mass, known 

The SED fitting in this work is mainly based on 
CIGALE FORTRAN where BC03 model is not included. BC03 
model is recently available in the CIGALE python version at 
http://cigale.lam.fr/pcigale/index_pcigale.html. We use this new 
version in testing the discrepancy caused by choice of synthesis 
models. 
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Figure 10. Top: SFR-M* relation of SPT DSFGs derived from CIGALE SED fitting using the single component SFH. The error bars are 
derived from the p robability distribution function by GIGALE. The dashed line show s the main sequence of star-forming galaxies at z = 2 
(IDaddi et al.ll20^b the d otted line is the relati on at z = 3 from IMagdis et HI (120101), and the dash-dotted line shows the relation at z =4 
at lower stellar masses bv IBouwens et ^ (120121 ). The medi an SFR — relat i on offGhapman et al.l (120051) SMGs at 2 ; = 2-3 is shown by 
the open diamond. These SMGs are extensively explored in IMichalowski et al.l (|2012l) by assuming different star formation histories which 
are reflected in the error bars. Bottom: The comparison between the stellar masses and SFRs derived using a single-component SFH and 
a double-component SFH. 


Table 4 

Input parameters for SED fitting with GIGALE. 


Parameter 

Symbol 

Values 

Metallicity 

Z 

0.02 (solar metallicity) 

Age in Gyr 

t 

0.01 0.05 0.1 0.5 1.0 2.0 

e-folding time r in Gyr 

T 

0.1 0.5 1.0 5.0 10.0 

V-band attenuation 

Av 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 

Slope correction of the Galzetti attenuation law 

5 

-0.2 -0.1 0.0 0.1 0.2 0.5 

IR power-law slope of dust mass distribution over heating intensity 

a 

1.2; 1.4; 1.5; 1.6; 1.7; 1.9 


as the “main sequence” (MS) of star-forming galaxies. 
The six SPT DSFGs presented in this paper all show 
intense star formation spanning a redshift range from 
2 ; = 2.8 to z = 5.7. We compare our results (Fig. [TQ|) 
wi th previous deterrn inations of the SFR — re l ations 
of iDaddi et al.l (|2QQ7D at 2 ; = 2 and iMagdis et al.l (|2Q1Q[ ) 
at 2 ; = 3 to investigate whether SPT DSFGs lie on top 
of the main-sequence or are extreme outliers. At 2 ; = 
2, IDaddi et al.l (|2QQ7[ ) derived the relation based on a 


sample of massive galaxies selected with K < 22 from 
GOODS and 24/im observations, w hich is widely used as 
the fiducial main sequence at 2 ; ^ 2. iMagdis et al.l (|2QlQf ) 
based their study on LBGs with IRAG observations, and 
found an increasing SFR at a given mass from 2 ; = 2 to 
2 ; = 3. All the stellar masses and SFRs are calibrated to 
a Ghabrier IMF (similar to Kroupa) and derived using 
the BC03 models. We do not correct our SMGs for the 
small differences between BC03 and M05 in the figure. 
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Figure 11. The SED fitting results for the six SPT DSFGs using the SED fitting code CIGALE, based on the assumptions of the M05 stellar 
population synthesis models and a single-component SFH. The data points from left to right are HST/WFG3 F110W+F160W/F140W, 
Spitzer/IRAG 3.5/rm+4.6/rm, Herschel/PAGS lOO^m+160/rm, Herschel/SPIRE 250^m+350/im+500/rm, APEX/LABOGA 87/rm, and SPT 
1.4mm+2.0mm. The circles denote detections while the triangles are 3(7 upper limits. 
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Figure 12. An intrinsic average SED (black) obtained by taking 
the median value of the six best-fit SEDs (color-coded) at each 
rest-frame wavelength. They are normalized to the total infrared 
luminosity of the average SED. 

SPT0538-50, SPT0103-45, and SPT2134-50 at ^ - 3 
lie above the z = 3 main sequence by factors of 3, 10, 
and at least 3, respectiv e ly. In contrast, the extensively 
studied I Chapman et alJ (j2QQ5f ) SMGs at z = 2-3, even 
though known to have widely varying stellar masses de¬ 
pending upon different assumed star formation histories, 
suggest that most SMGs represent the high-mass exten- 
si on of the main-sequence o f star-forming galaxies at 2; > 
2 (|Michalowski et al]l2Q12f ). Here we show the results as¬ 
suming the M05 models for direct comparison with our 
SMGs. Gomparing the median values of the two sam¬ 
ples, SPT SMGs have a factor of ^ 2 higher SFRs and a 
factor of ^ 2 lower stellar masses than those of the Ghap- 
man05 SMGs. This discrepancy still holds in the case of 
a double-SFH model (Fig. [TQ|) . This is suggestive of 
heterogeneity between SMG populations or two star for¬ 
mation modes: more quiescently star-forming SMGs and 
starburst SMGs (more discussion later in this section). 
We also show the individua l SMGs from the ALESS sam¬ 
ple (|da Cunha et al.ll2Q15f ) which is not a uniform pop¬ 
ulation: a significant fraction of the SMGs he above the 
MS while some are consistent with being at the high-mass 
end of the MS. 

6.2. Specific SFRs 

The specific SFR {SFR/Mfi)^ which can be interpreted 
as the inverse of the characteristic timescale for star for¬ 
mation to build the galaxy, is now recognized as a key 
parameter linked with the evolutionary status of these 
galaxies. Another way to explore the SFR — rela¬ 
tion is to look at sSFR as a function of redshift. The 
sSFRs for the six sources are presented along with stel¬ 
lar mass and SFR estimates in Table [5l We note that 
sSFRs are independent of lensing magnifications as long 
as differential magnification is not a concern. 

We compare our results wi th literature on th e sSFR- 
— z relations (Fig. [131 see iHeinis et al.l l2Q14l and ref¬ 
erence s therein). We show the results of iKarim et al.l 
jSonl), which span a wide s tellar mass range a t 0 < 2: 
< 3, and measurements from IHeinis et al ] (Iml . which 
extend from 2: = 1.5 to 4. These studies establish an 


upper boundary on the sSFR-z relations of normal star¬ 
forming galaxies in the literature. The dashed line is 
a fitting formula of the main-sequence sample of star- 
forming galaxies up to 2: = 4 in the COSMOS field by 
iBethermin et al.l (j2Q15a[ ) . The relations from these stud¬ 
ies enable us to make a fair comparison given that we 
make the same assumptions for IMF, SFH, and dust at¬ 
tenuation law in deriving the stellar masses, except that 
they adopt the BC03 model while we use M05. We stress 
that this ha s negligible e f fect o n our conclusions. 

Following IHeinis et al.l (I2Q14D. we demonst r ate th e ex- 
trapolation of the results of iBonwens et al.l (j2Q12f ) and 
iStark et~al] (j2Q13f ) from 2: = 4 to 7 at the relevant stel¬ 
lar mass bins by assuming that a power-law relation 
holds between SFR and at z > A but with a vari¬ 
ation in the slope between 0.7 and 1. We have seen 
the anti-correlation between sSFR and for the star¬ 
forming samples. SPT0103-45 and SPT2147-50 have sig¬ 
nificantly higher sSFRs, more than a factor of 10 above 
the MS, satisfying the c riterion for strong starbn rsts of 
>10 sSFRms defined in lBetherniin et aTI (j2015a|) . They 
significantly deviate from the main-sequence star forma¬ 
tion mode. SPT0538-50 lies about a factor of 4 above the 
MS and SPT2146-55 has a sSFR about 4 times higher 
with the lower limit reaching the MS. For SPT0346-52 
and SPT2134-50, we only have lower limits that already 
exclude MS SFRs. The mean sSFR for the four detec¬ 
tions lies above the MS at ^ 5cr confidence level. SPT 
SMGs, on average, ha ve sSFRs ^ 4 times h igher than 
those of ALESS SMGs (jda Cunha et al.l[2015f ) although a 
significant fraction of them are high-sSER outliers. These 
galaxies are not a uniform population. 

Hydrodynamical sim ulations with dust rad iative trans¬ 
fer calculations (e.g., iHavward et al.l l2Q12f ) have sug¬ 
gested the ‘bimodality’ of the SMG population: a mix 
of quiescently star-forming and starburst galaxies. The 
quiescent SMGs are star-forming disks blended into one 
submillimeter source at early stages of a merger whereas 
at late stages (near coalescence) tidal torques drive gas 
inflows and trigger strong starbnrsts. We are likely wit¬ 
nessing the ongoing strong starburst events. 

So far we have been applying the same lensing mag¬ 
nification factors derived from ALMA (i.e., dust emis¬ 
sion) to stellar emission. However, the effect of dif¬ 
ferential lensing between different wavelength regimes 
can be prominent in galaxy-galaxy lensing due to spa¬ 
tial variations within the background galaxy. Simula¬ 
tions ha ve predicted this effect in galaxy-g alaxy SMG 
systems (jHezaveh et al.l I2Q12I: iSerjeantI l2Q12f ) but we do 
not have the c apabi lity to measure it for our sample. 
iCalanog et al.l (j2Q14f ). which presents the NIR-derived 
lens models of Hcrschcl-^eXected galaxies, tests for dif¬ 
ferential lensing between the stellar and dust compo¬ 
nents. They find that the dust magnification factor mea¬ 
sured at 880 jam is ^ 1.5 times o n average higher than 
the near-IR magnif ication factor. iDve et all (j2Q15f ) and 
iRvbak et al.l (j2Q15[ ) find a similar differential magnifica¬ 
tion factor, based on the analysis of the ALMA high res¬ 
olution imaging of SDP.81. Assuming our SPT DSEGs 
have similar differential lensing effect, the resultant sS- 
ERs would decrease by ^1/3 which are within the error 
bars. 

Models predict that the SPT selection should identify a 
diverse set of sources. However. IBethermin et all (j2Q15bl ) 
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predicts that 90% of SPT leased sources are classified as 
main sequence galaxies. Even though the sample used 
for this study is small, the presence of extreme starbursts 
is thus surprising. One contributing factor may be the 
magnifications of the sources included in this subsample. 
Four of our sources have /i < 7 , while the average /i of t he 
full S PT sample predic t ed by IB ether min et al.l (j2Q12f ) is 
^15. iBethermiu et al.l (|2Q12f ) predict that sources with 
4 < /i < 7 on average lie a factor of 4.5 above the main 
sequence, compared to a factor of 1.7 for sources with 

/i > 20. 

7. SUMMARY AND CONCLUSIONS 

We have determined stellar masses and star forma¬ 
tion rates for six strongly lensed DSFGs using ALMA, 
HST^ and Spitzer follow-up observations. The sources 
are drawn from the SPT DSFG program. 

We develop a de-blending technique to extract emission 
of the background DSFG when blended with that from 
the foreground lensing galaxy in Spitzer/IRAC imaging. 
HSTWFC3 images are used to model the light profiles of 
the lenses which are subtracted in IRAC to produce the 
residual map overlaid with ALMA submillimeter imag¬ 
ing, indicating the position and structure of the gravita¬ 
tionally lensed DSFG. We demonstrate that in a gravita¬ 
tional lensing system like SPT0538-50, where the back¬ 
ground DSFG is strongly magnified and distorted by a 
massive lensing halo to well separate lensing arcs as seen 
by ALMA, a full Einstein ring structure is recovered in 
the residual map in both IRAC bands. Sufficient angu¬ 
lar separation is critical in terms of successfully disen¬ 
tangling the mixed emissions from the lens and source. 
This is best illustrated in SPT2134-50. In this case, there 
exists flux coincident with the ALMA contours, but we 
are unable to exclude the possibility that this signal is 
simply residual emission associated with the foreground 
lens. A minimum separation, comparable to the FWHM 
of the PSF, is required to conduct the de-blending and 
extract robust fluxes. 

Based on the multi-wavelength photometry, we have 
performed SED fitting with GIG ALE. This fitting in¬ 
cludes stellar synthesis models as well as dust emission 
templates (i.e., taking the full spectral energy distribu¬ 
tion into consideration) to derive stellar masses, total 
IR luminosities, and instantaneous star formation rates. 
The best-fit model is based on a single extended burst 
SFH. The derived stellar masses span a wide range with 
a median value of ^ 5 x IO^^Mq. The instantaneous 
SFRs range from 510 to 4800 M0yr“^. Use of a two- 
component SFH generally leads to higher stellar masses, 
higher SFRs but similar sSFRs, and our qualitative re¬ 
sults are robust to choice of SFHs. 

We investigated the SFR — relation and evolu¬ 
tion of specific SFR with redshift. Compared to sSFR 
- z relations of normal star-forming galaxies in the lit¬ 
erature, the six SPT DSFGs all lie above the MS. Two 
of them are ^ 10 times above the star-forming main- 
sequence. SPT0346-52 converts the gas into stars at a 
rate of ^ 4800 Mqyt~^ only 1 Gyr after the Big Bang, 
a SFR that is among the highest at any epoch. Our 
results suggest that we may be witnessing strong star- 
burst events. One possible explanation for such extreme 


starbursts is gas-rich major mergers (e.g 

iTacconi et al. 

120061: iNaravanan et al.lDOOOUEngel et al.l 

l2010l:lDve et al. 


I2015D in which tidal torque compresses gas reservoir lead¬ 
ing to a concentrated starburst. High resolution kine- 
matical and morphological observations are needed to 
test the merger scenario. 

Future work involves determining dynamical mass 
from CO line width and gas mass from CO luminosity 
together with dust mass estimates to provide a more 
comprehensive interpretation of this unique population 
(Spilker et al. a nd Strandet et al. in prepara tion). As 
noted by others (jMichalowski et al1l2Q 1^1201^ . deriving 
stellar masses for young, dusty DSFGs is an exercise 
fraught with challenge. For high redshift sources, 
the lack of rest-frame NIR bands further limits the 
achievable fidelity. Observations with JWST offer the 
potential for dramatic improvement due to the two-fold 
advantage of improved spatial resolution and greater 
wavelength coverage extending into the mid-infrared. In 
the MIR, which probes rest-frame Ks for distant SMGs, 
the foreground lens fades while the SMG brightens. The 
improved spatial resolution compared to Spitzer greatly 
simplifies separation of source and lens fluxes, even for 
compact lensing configuration. Moreover, the required 
integration times will be short. For instance all the 
sources are expected to be detected within 1 minute 
integration time at S/N > 10 at rest-frame Ks band. 

ACKNOWLEDGMENTS 

We thank the anonymous referee for useful comments. 
We thank Sebast ien Heinis for sharing the data in 
iHeinis et al.l (j2Q14f ). This work is based on observations 
made with the NASA/ESA Hubble Space Telescope, ob¬ 
tained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with programs 7^12659 
and #13614. Support for programs #12659 and #13614 
was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc., 
under NASA contract NAS 5-26555. This work is based 
in part on observations made with the Spitzer Space Tele¬ 
scope, which is operated by the Jet Propulsion Labo¬ 
ratory, California Institute of Technology under a con¬ 
tract with NASA. Support for this work was provided by 
NASA through an award issued by JPL/Caltech. Par¬ 
tial support for this work was provided by NASA from 
awards for Herschel observations for OTl_dmarrone_l, 
OTl_jvieira_4, and OT2_jvieira_5. Herschel is an ESA 
space observatory with science instruments provided by 
European-led Principal Investigator consortia and with 
important participation from NASA. J.S.S., D.P.M., and 
J.D.V. acknowledge support from the U.S. National Sci¬ 
ence Foundation under grant No. AST-1312950 and 
through award SOSPAl-006 from the NRAO. M. A. ac¬ 
knowledges partial support from FONDECYT through 
grant 1140099. This paper makes use of the following 
ALMA data: ADS/JAO.ALMA #2011.0.00957.S and 
#2011.0.00958.S. ALMA is a partnership of ESO (rep¬ 
resenting its member states), NSF (USA) and NINS 
(Japan), together with NRG (Canada) and NSC and 
ASIAA (Taiwan), in cooperation with the Republic of 
Chile. The Joint ALMA Observatory is operated by 




















17 


Table 5 

Derived properties from the CIGALE SED fitting 


SPT source name 

stellar mass 

L/fl [8 - 1000 

age 

SFRcigale 

SFRir 

sSFR 

ID 

Mq 

L© 

Myr 

Mo/yr 

Me/yr 

Gyr-i 

SPT0103-45 

5.515:9 ^ 10^° 

(1.2 ±0.1) X lO^® 

40155 

17401555 

1290-1800 

SI 4 + 40.1 

SPT0346-52 

< 3.1 X 10“ 

(3.6 ±0.3) X 10i3 

3015° 

484011°°° 

3830-5340 

> 15.7 

SPT0538-50 

3.9l5;g X 10^° 

(4.3 ±0.5) X IOI 2 

or »-|-250 

SlOllfo 

460-640 

13.311+ 

SPT2134-50 

< 6.6 X 10i° 

(4.9 ±0.7) X IOI 2 

40+"^® 

640l“° 

520-730 

> 9.1 

SPT2146-55 

0.8lJ;g X 10 “ 

(9.3 ±1.2) X 10^2 

70+250 

1190155° 

990-1380 

.. 0+54.0 
lO.O_i2 2 

SPT2147-50 

2.0lJ;| X 10i° 

(8.4 ±1.0) X 10^2 

4U io 

1290152° 

900-1250 

64.ol°|;° 


Note. — The derived properties are corrected for leasing magnifications. For SPT2134-50 and SPT0346-52 which are non-detections in 
the IRAC bands, we place 3(7 upper limits on the stellar masses. We note that specific SFRs are independent of magnification factors. 


10"^ ^ ^ ^ ^ ^ ■ r 

★ SPTSMGs 

ALESS SMGs (da Cuhna+2015) 

• Negrello+2014 


single-component SFH 



Karim+2011 logM* = 9.8-10.2 
Karim-H2011 logM* = 10.2-10.6 
Karim-K2011 logM* = 10.6-11.0 
Karim-H2011 logM* >11.0 
Bouwens-i-2012 logM* > 10.0 
Stark-i-2013 logM* > 10.0 
Heinis-i-2014 logM* = 10.0 10.5 
Bethermin-i-2015 logM* > 10.4 


11.0 = 


Figure 13. sSFR - z relation. The red filled stars are the SPT DSFGs. SPT2134-50 and SPT0346-52 have lower limits on sSFRs 
(SPT2134-50 is shifted by -0.08 in redshift for clarity). We show the results from IKarim et a ll (|20^ ) (green ban ds) which cover a board 
range of redshift and stellar mass bins. The blue shaded regions are average sSFRs obtained bv [ Heinis et al.| ([2014|1 at z 1.5, 2 ; ~ 3, and 2 ; 
~ 4. The pink and orange shaded regions are the extrapolated relations bv IBouwens et af ] (I 2 OI 2 I) and lStark et al.l (|2 013l). The black curve 
is a fitting formula of the main-sequence sample of star-forming galaxies up to 2 ; = 4 in the GOSMOS field bv IBethermin et al.l (I2015bli . The 
magenta squares are individua l lensed DSFGs from the literature with sSFRs high abo ve the MS: Gos mic Eyelash (i.e., SMMJ 2135-0102) at 
z = 2.3 (ISwinbank et al.ll2010l ). HXMMOl at z = 2.3 (shifted by -0.15 in 2 ; for clarity: TFu et al.ll2013l ). HLSW-01 at 2 ; = 2.96 (IGonlev et al.l 
[20Tlh . HATLAS12-00 at z = 3.2 (IFn et al.ll20Tli . AzTEG-3 at z = 5.3 (IGaoak et al.ll201in. and HF LS3 at z = 6.34 (IGoorav et al.ll2014l) . 
The magenta circles are SDP.17, SDP.81, and SDP.130 from Herschel-ATLAS (jNegrello et al.ll2014l i. The grey dots are 2 : > 2 SMGs from 
the ALESS sample (jda Gunha et al.ll201^ b SPT0103-45 and SPT2147-50 have si gnificantly higher sSFRs, more than a factor of 10 above 
the MS, satisfying the criterion for strong starbursts of > 10 sSFRms defined in IBethermin et al ] (I2015bl b 
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